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ABSTRACT 



> 



We present Chandra and XMM-Newton X-ray observations that monitor the neu- 
tron star cooling of the quasi-persistent neutron star X-ray transients KS 1731—260 
and MXB 1659—29 for approximately 4 years after these sources returned to quies- 

■ cence from prolonged outbursts. In both sources the outbursts were long enough to 
| significantly heat the neutron star crust out of thermal equilibrium with the core. 
, We analyse the X-ray spectra by fitting absorbed neutron star atmosphere models to 

■ the observations. The results of our analysis strengthen the preliminary findings of 
| Wijnands et al. that in both sources the neutron star crust cools down very rapidly 
. suggesting it has a high heat conductivity and that the neutron star core requires en- 
hanced core cooling processes. Importantly, we now detect the flattening of the cooling 

^ (— | in both sources as the crust returns to thermal equilibrium with the core. We measure 

O ^ the thermal equilbrium flux and temperature in both sources by fitting a curve that 

decays exponentially to a constant level. The cooling curves cannot be fit with just a 
simple exponential decay without the constant offset. We find the constant bolometric 
flux and effective temperature components to be (9.2 ±0.9) x 10~ 14 ergs cm -2 s _1 and 
70.0 ± 1.6 eV in KS 1731-260 and (1.7 ±0.3) x 10~ 14 ergs cm" 2 s" 1 and 51.6 ± 1.4 eV 
in MXB 1659—29. We note that these values are dependent on the assumed distance 
to the sources and the column density which was tied between the observations due to 
the low number of photons in the latter observations. However, importantly, the shape 
j^j ■ of the cooling curves is independent of the distance assumed. In addition, we find that 

the crust of KS 1731-260 cools faster than that of MXB 1659-29 by a factor of -2, 
likely due to different crustal properties. This is the first time that the cooling of a 
neutron star crust into thermal equilibrium with the core has been observed in such 
detail. 

Key words: accretion, accretion discs — X-rays: binaries — stars: neutron — X-rays: 
individual (KS 1731-260), X-rays: individual (MXB 1659-29) 
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1 INTRODUCTION 

Low-mass X-ray binaries consist of either a neutron star or 
a black hole that accretes matter from a low-mass compan- 
ion star, where mass transfer occurs because the donor star 
overflows its Roche-lobe. The soft X-ray transients form a 
sub-group of low-mass X-ray binaries. The vast majority 
of their time is spent in a quiescent state where little ac- 
cretion occurs, resulting in X-ray luminosities < 10 34 ergs 
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s _1 . However, during outbursts (which usually last weeks 
to months) the mass-accretion rate rises significantly, with 
a corresponding large increase in brightness - the typical 
outburst luminosity is in the range 

10 36 _ in 39 erg s" 1 (for 
reviews of neutron star and black hole transients see , e.g., 
ICampana et alJll998allRemillard fc McClintocldl2006l) . 

The quiescent X-ray spectra of neutron-star X-ray tran- 
sients can generally be characterised by two components 
- a soft, thermal component that dominates below a few 
keV, and a hard power-law tail that dominates above a 
few keV. The power-law component is not always present 



© 0000 RAS 



2 Cackett et al. 



and its origin is poorly understood, though in some sys- 
tems it completely dominat es (e.g. EXO 1745—248 in the 
globular cluster Terzan 5, IWiinands et alJ l2005al) . Sev- 
eral emission mechanisms have been proposed to explain 
the quiescent X-ray emission, for example, residual accre- 
tion onto the neutron star magnetosp h ere or pulsar shock 
emission (e.g.. | C ampana et al.lll998at iMenou et all Il999t 
ICampana fc Stelrall200ri iMenou fc McClintockll200ll) . How- 
ever, the mechanism most often used to explain the thermal 
X-ray emission is that in which the emission is due to the 
cooling of a hot neutron star which has been heated d uring 
out burst ivan P aradiis et aljll9 8^: iBrown et alJll99§t) . In 
the lBrown. Bildstoi^^Rutiedge) (I1998T) model, the neutron 
star core is heated by nuclear reactions occurring deep in the 
crust due to accretion during outburst, and this heat is re- 
leased as thermal emission during quiescence. Thus, the tem- 
perature of the neutron star core and surface in quiescence 
is then set by the time-averaged mass accretion rate and 
the neutron star properties. Therefore, if the time-averaged 
accretion rate can be estimated, studying quiescent neutron 
star systems can give important information about the phys- 
ical processes at work in the very high densities present in 
neutron star cores. 

Some X-ray transients (the quasi-persistent transients) 
spend an unusually long time in outburst, with outbursts 
lasting years to decades rather than the more common 
weeks to months. For the neutron-star quasi-persistent tran- 
sients, during outbursts, the crust of the neutron star is 
heated to a point beyond thermal equilibrium with the 
stellar interior. Once accretion falls to quiescent levels, the 
crust cools thermally, emitting X-ra ys, until it reaches equ i- 
librium again with the core (e.g., iRutledee etUl Eo02b ) . 
Crust cooling curves allow us to probe the properties of 
neutron stars, such as, the kind of core cooling processes 
at work (which depends on the equation of state of ul- 
tra dense matter) and the heat conductivity of the crust 
(which likely depends on the properties of the iron lat- 
tice of the crust, such as its purity). While it is ex- 
pected that the crust can cool significantly between out- 
bursts, the neutron star core should only cool significantly 
on much longer timescales. X-ray observations of the two 
quasi-persistent neutron s tar X-ray transie n ts KS 1731—260 
JWiinands et alJ l2001al: IWiinandsl l2002t IWiinands et alJ 



l2002al : iRutledee et alJ l2002bj : IWiinandsM2005l) and MXB 
1659-29 iWiinands et alJl2003l . l2004l : IWiinandsll2004 120051) 

have tracked the initial cooling of these neutron stars once 
they went into a quiescent state. In this paper, we present 
additional Chandra and XMM-Newton observations that 
further track the cooling of these two neutron stars as they 
return to thermal equilibrium with the core. 

1.1 KS 1731 260 

KS 1731—260 was first detected in outburst in August 1989 
by the Mir-Kvant instrument and further analysis revealed 
that t he source was also found to be act ive in October 
1988 iSunvaevlll989l: ISunvaev et all Il990ali . As no previ- 
ous telescopes had detected this source, it was classified 
as a soft X-ray transient and the detection of type I X- 
ray bursts identified the compact object as a neutron star. 
Since then, until early 2001, the source remained bright 
and was observed by many X-ray telescopes including the 
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Figure 1. RXTE All-Sky Monitor 7-day averaged lightcurve for 
KS 1731-260 (top) and MXB 1659-29 (bottom). The times of 
the Chandra and XMM-Newton observations of these sources 
during quiescence are marked by solid and dotted vertical lines 
respectively. The last detection of KS 1731—260 in outburst 
was 1847 days after 1996 January 1, and the last detection of 
MXB 1659-29 was 2076 days after 1996 January 1. 
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In early 2001, KS 1731 — 260 was undetected by monitoring 
observations by the All-Sky Monitor (ASM) on-board the 
Rossi X-ray Timing Explorer (RXTE) as can be seen by 
the ASM lightcurve of this source in Fig. (top). Further 
pointed RXTE observations with the Proportional Counter 
Array (PCA) and the bulge scan observations of the Galac- 
tic centre region did not detect the source, with the first 
non-detection by the PC A observations on 2001 February 7 
JWiinands et alJl2001al) . Using archival RXTE/PCA data, 
we find that the final detection of the source with RXTE 
was on 2001 January 21, thus it was still in outburst and 
actively accreting on this date. We extracted the RXTE X- 
ray spectrum of the source on that day and find that a 2-10 
keV unabsorbed flux of ~ 2.3 x 10 -10 ergs cm -2 s _1 which 
corresponds to a luminosity of ~ 1 X 10 36 erg s _1 for a 
distance of 7 kpc (we assume a distance to KS 17 31—260 
of 7 kpc throughout the paper, iMuno et al1l200d) . Thus, 
KS 1731 — 260 went into a quiescent state sometime after 
2001 January 21 and it was actively accreting for around 
12.5 years. This last detection in outburst gives a much bet- 
ter constraint than had previously been found for when this 
source turned off. Unfortunately, there are no other known 
outbursts of this source, thus we are unable to tell whether 
this is a typical outburst for this source, or whether previous 
outbursts have been significantly longer or shorter. 

Within two months of KS 1731-260 going into a qui- 
escent phase there was a ~20 ks Chandra Advanced CCD 
Imaging Spectrometer (using AC IS-S mode) observatio n 
of this source on 2001 March 27 llWiinands et al]l2001al) . 
These authors detected the source at an unabsorbed flux of 
~2x 10~ 13 ergs cm -2 s _1 (0.5-10 keV), which corresponds 
to a 0.5-10 keV luminosity of ~ 1 x 10 33 ergs s _1 , a factor 
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of 10 3 lower than at the last RXTE/PCA detection. Obser- 
yations by BeppoSA X confirmed the low flux of the source 
jBurderi et all (2002ft re ported by IWiinands et all fcOOlaft . 
IWiinands et al.l l)2001aft note that this luminosity is much 
lower than expected if this is due to the cooling of the neu- 
tron star unless the source spends an unusually long amount 
of time in a quiescent state (at least several hundred years) 
and even longer if the crust is dominating emission from 
the neutron star, or that alternatively, enhanced cooling of 
the core occurs due to higher levels of neutrino emission. 
IRutledee et al ] J2002bft note that given the long outburst, 
the crust would be dominating the thermal emission from 
the neutron star as it would have been heated significantly 
out of equilib rium with the core. Giv en the accretion history 
of this source IRutledee etafl J2002bft simulate the quiescent 
lightcurve for the source, finding that depending on the crust 
and core microphysics, the cooling timescale for the crust to 
return to equilibrium with the core ranges from 1 to 30 years. 

Additional X-ray observations of this source in quies- 
cence with XMM-Newton in 2001 September 13-14 revealed 
that the source had dimmed further to a 0.5- 10 keV luminos- 
ity of approximately (2 — 5) x 10 32 ergs s _1 llWiinands et alJ 
l2002aft . Under the neutron star cooling model, this decrease 
in luminosity is interpreted as rapid cooling of the the neu- 
tron star crust (a factor a ~2-3 in just half a year) suggesting 
that the crust has a high thermal conductivity and that en- 
hanced core cooling processes are required for the core to be 
so cool. Preliminary analysis of two extra Chandra observa- 
tions of this source on 2003 February 11 and 2003 August 8 
yet again showed a further decrease in the luminosity of the 
source with a possible flattening of the luminosity cooling 
curve suggesting that t he crust may hav e returned to equi- 
librium with the core JWiinanddl2005ft . In this paper we 
present the final analysis of these Chandra data, re-analyse 
all the previous Chandra and XMM-Newton observations of 
this source in quiescence and also present 2 additional Chan- 
dra observations and 1 additional XMM-Newton observation 
of this source. 

1.2 MXB 1659 29 

MXB 1659—29 is an X-ray transient discovered by 
iLewin et al] i ll 9 76ft with SAS-3 in 1976 October. Type-1 X- 
ray bursts detected by these initial observations established 
that the compact object was a neutron star. Between 1976 
October and 1978 September, this source was detected sev- 



eral times bv both SAS-3 and HE AO 1 jLewin et al 
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Griffiths et al. 19781: iDoxsev et al. 


197S 


Cominskv et alJll983HCominskv & Woodlll984. 1989 


). The 



source exhibits X-ray eclipses a nd dips with a period of 
~ 7.1 hr, due to the orbital period llCominskv fc Woodll984l 
1989). In 1 979 July, the Hakucho satellite failed to detect 
the source llCominskv et al] |l983ft, and later observations 
by EX OSA rand ROSAT llOosterbroek et alJl200ll : IVerbuntl 
2001) failed to detect the source also. The 0.5-10 keV upper 
limit on the unabsorbed flux from the ROSAT obs ervations 
in 1991 and 1992 is (1 - 2) x 10~ 14 ergs cm" 2 s" 1 llVerbuntl 
1200 ll: IWiinandsll200a lOosterbroek et aljElOlft . 

The source remained in a quiescent state for a period 
of 21 years, until on 1999 April 2 BeppoSAX d etected the 
source in outburst again (|in 't Zand et al]ll999ft . This out- 
burst was also detected by the RXTE ASM (see the ASM 



lightcurve in Fig. and studied in detail by pointed obser- 
vations with the RXTE PCA, BeppoSAX and XMM-Newton 
llWachter et al] l200d lOosterbroek et al.lEoOll: ISidoli et al] 
l200lt IWiinands et aLlEoOlbl l2002cft~ 

After an outburst of approximately 2.5 years in length, 
MXB 1659—29 went back into a quiescent state in 2001 
September. The last detection of the source was on 2001 
September 7 with the RXTE PCA, and the source was 
not de tected on 2001 September 14, 24, 30 iWiinands et al] 
l2002cft . Just over a month after the source returned to 
a quiescent state, Chandra observations, on 2001 Octo- 
ber 15, detected it with an unabsorbed 0.5-10 keV flux of 
(2.7 — 3.6) x 10~ 13 ergs cm -2 s^ 1 which corresponds to a 
0.5-10 keV lumino sity of (3.2 - 4.3) x 10 33 ergs s" 1 for a 
distance of 10 kpc llWiinands et al.|l2003ft . This flux is a fac- 
tor greater than 10 higher than the ROSAT non-detection 
upper limit during quiescence in the early 1990s, suggesting 
that the emission is dominated by thermal emission from 
the neutron star crust that has been heated significantly 
out of thermal equilibrium with the core. However, as in 
the case of KS 1731—260, the flux was lower than expected 
given the accretion history, suggesting that enhanced cool- 
ing processes are occurring in the core. Further Chandra 
observations of MXB 1659—29 in quiescence on 2002 Octo- 
ber 15 and 2003 May 9 showed that the flux was decaying 
exponentially, with the source flux decreasing by a factor of 
7-9 over the 1.5 years between the 2001 and 2003 Chandr a 
observations llWiinands et al]l2004 Iwiinandsll200i l200Sft . 
This is consistent with the accretion-heated crust cooling 
down in quiescence. The rapid timescale of the cooling im- 
plies that the crust may have a large therm al conductivity 
llWiinands et aJ]l2004lWiinanddl20olll200Sft . In this paper 
we re-analyse all the previous observations of this source in 
quiescence with Chandra and also present 2 additional Chan- 
dra observations and 1 additional XMM-Newton observation 
of this source. 



2 OBSERVATIONS AND ANALYSIS 
2.1 KS 1731 260 

Here we analyse 5 Chandra and 3 XMM-Newton observa- 
tions of KS 1731 — 260 in a quiescent state spread over ~4 
years since the end of the outburst in 2001 February. Details 
of all the observations are given in Table We will first 
describe the analysis of the Chandra and then the XMM- 
Newton data. 



2.1.1 Chandra analysis 

All the Chandra observations were taken in ACIS-S con- 
figuration. The Chandra data were reprocessed and anal- 
ysed using CIAO (v3.3) and Caldb (v3.2.1) following the 
standard analysis threads. For all of the Chandra obser- 
vations, the source lightcurve and spectrum was extracted 
from a circle of radius 3" around the source position, and 
the background lightcurve and spectrum was extracted from 
a source-free annulus with inner radius 7" and outer radius 
25". We checked the background lightcurve for significant 
background flares, and none were found. The background- 
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Table 1. Details of the Chandra (CXO) and XMM-Newton (XMM) observations of KS 1731-260. The XMM-Newton observations 
0137950201 and 0137950301 were so close in time that they were merged together. Here we give the details for the merged data. The 
background-subtracted net count rate is for the 0.5-10 keV band. 



ObsID 


Date 


Good Time 


Net Count Rate 


(Telescope) 




(ks) 


(counts s — 1 ) 




9001 l\/Tcirr>Vi 97 


1 Q A 

_i_y .44 


fl + (1 7W 1 n— 3 


0137950201/301 (XMM) 


2001 September 13/14 


46.0 (MOS1) 


(1.7 ±0.2) x 10~ 3 (MOS1) 






46.1 (MOS2) 


(2.0 ±0.2) x 10- 3 (MOS2) 






33.8 (pn) 


(3.7 ±0.4) x 10" 3 (pn) 


3796 (CXO) 


2003 February 11 


29.7 


(1.4 ±0.2) x 10" 3 


3797 (CXO) 


2003 August 8 


29.7 


(1.1 ±0.2) x 10~ 3 


0202680101 (XMM) 


2005 February 28/March 1 


70.0 (MOS1) 


(2.4 ±1.6) x 10" 4 (MOS1) 






70.0 (MOS2) 


(6.4 ±1.7) x 10" 4 (MOS2) 






46.0 (pn) 


(1.1 ± 0.4) x 10" 3 (pn) 


6279 (CXO) 


2005 May 10 


10.4 


(0.5 ±0.3) x 10~ 3 


5468 (CXO) 


2005 June 4 


38.5 


(0.8 ±0.2) x 10- 3 



subtracted net count rates in the 0.5-10 keV band are given 
in Tabled 

2.1.2 XMM-Newton analysis 

The first two XMM-Newton observations were separated 
by only a matter of hours, therefore we c ombine both of 
these t ogether to increase sensitivity (as did lWiinands et alJ 
2002a). Here we analyse data from the three EPIC instru- 
ments on-board XMM-Newton - the two metal oxide semi- 
conductor (MOS) cameras and the pn camera which were 
operated in full-frame mode with the thin filter in all the 
observations. To analyse the XMM-Newton data we use the 
Science Analysis Software (SAS, version 6.5.0). We extract 
source lightcurves and spectra from a circle of radius 10". 
The background extraction region was chosen to be a source- 
free circle of radius 1' close to the source position. An an- 
nulus was not used because of the location of a nearby 
source, 2MASSI J173412. 7-260548, approximately 30" away. 
Although this source is also present in the Chandra data the 
significantly better angular resolution of Chandra allows the 
use of an annulus for background extraction. 

In the merged observations (ObsID 0137950201/301) 
background flares were present. We kept all data where the 
count rate for PI > 10 keV was less than 10 counts s _1 for 
both the MOS and pn. In the last XMM-Newton observation 
(ObsID 0202680101) background flares were also present. In 
this case we kept all data where the count rate for PI > 
10 keV was less than 10 counts s _1 for the MOS and less 
than 15 counts s _1 for the pn. We note that in this last ob- 
servation the source is not obviously present when viewing 
an image in the full energy range. It only becomes appar- 
ent when restricting the energy range to between 1 and 2 
keV where the vast majority of the photons from the source 
are emitted during this observation. The net count rates for 
these observations are given in Table Q 

2.1.3 Spectral analysis 

For spe c tral analysis of the data we use xspec 
iArnaudl Il996l) . The data was fit with an ab- 
sorbed neutron star hydrogen atmosphere model 



(nonmagnetic case, IZajdin. Pavlov fc Shibanovl Il996l : 
iPavlov. Shibanov fc Zavlinlll99ll) with the mass and radius 
fixed at the canonical values (1.4 Mq and 10 km respec- 
tively). Both the Chandra and XMM-Newton data were fit 
simultaneously with the column density tied between the 
observations and the temperature left as a free parameter. 
The column density was tied between the observations 
as there were too few counts in the last observations to 
constrain this parameter. The parameters for the EPIC 
MOS1, MOS2 and pn spectra for each XMM-Newton 
observation were all tied. The normalisation of the neutron 
star atmosphere model was fixed for a di stance of 7 kpc 
(distance determined bv lMuno et alJ 12000). Although this 
distance is uncertain, if this parameter is left free it goes 
to an unrealistic value that would put the distance to the 
source at ~ 20 kpc (although the error on the distance is 
very large and includes the assumed 7 kpc distance). Hence, 
we chose to fix the distance at 7 kpc. 

The data were unbinned, due to the low number of 
counts in the last f ew observations, and the W statistic 
(Wa chter et al .119791) was used in the spectral fitting process 
with background subtr action. The W statistic is adapted 
from the Cash statistic <ICasblll979l) . but allows background 
subtraction. We note that spectral fitting to the first Chan- 
dra observation using both \ 2 an( l W-statistics resulted in 
the same parameters to within the errors. To determine the 
bolometric flux we extrapolated the model to the energy 
range 0.01 - 100 keV, which gives approximate bolometric 
fluxes. The results of the spectral fits are detailed in Ta- 
ble [5] The spectral parameters for the first Chandra and 
XMM-Ne wton observations are con sistent with those deter- 
mined by IWiinands et al ] J2002al) . We also note that we 
get results co nsistent with the previou sly published data for 
KS 1731-260 fWiinands et alJ j2002al) when using alterna- 
tive models to the NSA model (i.e., a blackbody model). 



2.1.4 Cooling curves 

It is clear from Table [5] that both the effective temperature 
for an observer at infinity (T^§) and the bolometric flux fboi 
decrease significantly with time. This cooling cannot be fit 
by a simple exponential decay, giving reduced % 2 values of 
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Table 2. Neutron star atmosphere model fits to the X-ray spectrum of KS 1731—260 for 5 Chandra (CXO) and 3 XMM-Newton 
(XMM) observations. Luminosity is calculated assuming a distance to the source of 7 kpc. 1-ct errors on the parameters are given. The 
Modified Julian Date (MJD) given correspond to the mid-point of the observation. Note the column density parameter is tied between 
the observations when simultaneously fitting the data, and we determine Nh = (1-3 ± 0.1) X 10 22 cm -2 . 



ObsID MJD kT^j Bolometric flux Luminosity 

(Telescope) (eV) (10~ 13 ergs cm -2 s" 1 ) (10 32 ergs s" 1 ) 



2428 (CXO) 


51995.1 


103^3 


4 O+0.3 


24+ 2 


0137950201/301 (XMM) 


52165.7 


88 ± 2 


2.2 ±0.2 


13 ±1 


3796 (CXO) 


52681.6 


76 ± 3 


1.2 ±0.2 


7± 1 


3797 (CXO) 


52859.5 


72 ± 4 


1.0 ±0.2 


6± 1 


0202680101 (XMM) 


53430.5 


70 ±3 


9 +0 - 2 


5± 1 


6279 (CXO) 


53500.4 


67tg 


0.8 ±0.3 


4±2 


5468 (CXO) 


53525.4 


70 ±3 


0.9 ±0.2 


5± 1 



4.3 and 4.8 for fits to the temperature and flux curves (see 
the dotted curves in Fig. [5] However, it is fit well by an 
exponential decay that levels off to a constant offset of the 
form y(t) = aexp [— (t — to)/b] + c, with a a normalisation 
constant, b the e- folding time, c a constant offset, and to the 
start time. When fitting to the data, to was fixed to midday 
on the last day that the source was observed to be active, 
MJD 51930.5, though we find that the other parameters are 
not very sensitive to the exact value of to- The best-fitting 
cooling curves are shown in Fig. H For the T£§ curve a = 
39.5±3.6 eV, b = 325±101 days, and c = 70.0±1.6 eV, with 
xl = 0.2. For the F bo i curve a = (4.3±0.5) x 10" 13 ergs cm" 2 
s" 1 , b = 212±60 days, and c = (9.2±0.9) x 10" 14 ergs cm" 2 
s" 1 , withxl = 0.4. The reason for the difference between the 
e-folding times of the flux and effective temperature curves 
is explained in Appendix 1X1 



2.2 MXB 1659 29 

We analyse 5 Chandra observations and 1 XMM-Newton ob- 
servation of MXB 1659—29 whilst the source was in a qui- 
escent state spanning a period of ~4 years after the end of 
the outburst in September 2001. Details of the observations 
are given in Table [3] We will first describe the Chandra and 
then the XMM-Newton data reduction and analysis. 



2.2.1 Chandra analysis 

All the Chandra observations of this source were taken 
in the ACfS-S configuration. As in the analysis of the 
Chandra data for KS 1731 — 260, we reprocess and anal- 
yse the data using CIAO (v3.3), Caldb (v3.2.1) and the 
standard analysis threads. For all of the Chandra obser- 
vations, the source lightcurve and spectrum was extracted 
from a circle of radius 3" around the source position, and 
the background lightcurve and spectrum was extracted from 
a source-free annulus with inner radius 7" and outer radius 
22". We checked the background lightcurve for significant 
background flares, and none were found. 

The analysis of the data for MXB 1659—29 is compli- 
cated by the fact that this source is ecli psing with an eclipse 
duration of ^900 s and period of 7.1 hr jC ominsky fc Woodl 
11984. 119891 : IWachter et aflfeOOCt lOosterbroek et al.l 120011) 
and so we receive no (or minimal at mo st) counts from 
the source during the eclipse in quiescence iWiinands et alJ 
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Figure 2. Cooling curves for KS1731— 260. Top: Bolometric flux 
versus time since to; the last detection of the source accreting. 
The solid line shows the best fitting exponential decay to a con- 
stant. The constant is indicated as a dashed line. The dotted line 
shows the best fitting simple exponential decay, which does not 
fit the data well. Bottom: Effective temperature for an observer 
at infinity versus time since the end of the outburst. The lines are 
as above. 



2003). While there are enough counts in the first Chan- 
dra observation to d etect the eclipse in the lightcurve 
JWiinands et alJl2003l) . this is impossible with later obser- 
vations and so we manually reduce the exposure time by 
900 s to compensate for the source being in eclipse, having 
checked that only one eclip se occurs during e a ch ob serva- 
tion using the ephemeris of lOosterbroek et alJ (J2001) . The 
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Table 3. Details of the Chandra (CXO) and XMM-Newton (XMM) observations of MXB 1659-29. The background-subtracted net 
count rate is for the 0.5-10 keV band. 



ObsID 


Date 


Good Time 


Net Count Rate 


(Telescope) 




(ks) 


(counts s L ) 


2688 (CXO) 


2001 October 15 


17.9 


(5.2 ±0.2) X 10" 2 


3794 (CXO) 


2002 October 15 


26.1 


(9.8 ±0.6) x 10" 3 


0153190101 (XMM) 


2003 March 13 


65.3 (MOS1) 


(3.2 ±0.3) x 10~ 3 (MOS1) 






65.4 (MOS2) 


(2.9 ±0.3) x 10" 3 (MOS2) 






47.1 (pn) 


(1.1 ± 0.1) x 10- 2 (pn) 


3795 (CXO) 


2003 May 9 


26.2 


(3.9 ±0.4) x 10" 3 


5469 (CXO) 


2005 July 8 


27.7 


(1.1 ±0.2) x 10~ 3 


6337 (CXO) 


2005 July 25 


17.8 


(0.7 ±0.2) x 10~ 3 



background-subtracted net count rates in the 0.5-10 keV 
band are given in Table [3] 

2.2.2 XMM-Newton analysis 

To analyse the XMM-Newton observation of MXB 1659-29 
we use SAS (version 6.5.0). We analyse data from the EPIC 
MOS1, MOS2 and pn cameras which were operated in full- 
frame mode with the medium filter (the source is too faint 
to be detected on the RGS). The source lightcurve and spec- 
trum was extracted from a circle with a 15"radius. The 
background lightcurve and spectrum was extracted from a 
source- free annulus with inner radius of 40"and outer radius 
of 100". Several background flares were present during the 
observation. We chose to keep all data where the count rate 
for PI > 10 keV was less than 7 counts s" 1 for the MOS 
and less tha n 10 counts s~ for the p n detector. Using the 
ephemeris of lOosterbroek et alJ feOOll) we determined that 3 
eclipses occurred during the good time of the MOS observa- 
tions, and so we manually reduce the exposure times for the 
MOS observations by 2.7 ks. In the pn observations several 
flaring events occurred during the eclipse times. We deter- 
mine the total amount of good time for the pn during the 
eclipses to be 1.675 ks, and thus we reduce the exposure time 
by this amount. The net count rates for this XMM-Newton 
observation are given in Table |3] 

2.2.3 Spectral analysis 

We perform a similar spectral analysis for the MXB 1659—29 
data as for the KS 1731-260. The Chandra and XMM- 
Newton spectra are left unbinned and fit simultaneously 
with an absorbed non-magnetic neutron star hydrogen at- 
mosphere model using the W statistic with background sub- 
traction. The column density is tied between the separate 
observations, the mass and radius fixed at their canonical 
values, and the temperature was the only parameter left 
free between the observations. The parameters are tied be- 
tween the MOS1, MOS2 and pn spectra of the XMM-Newton 
observat ion. The distance to this source is in the ra nge 10- 
13 kpc jMuno et alJl200ll: lOosterbroek et al.ll200ll) . When 
leaving the normalisation as a free parameter in the Chan- 
dra fits this leads to a distance to the source of 7.8 kpc with 
no constraint on the upper limit of the distance and a lower 
limit of 6.2 kpc. So that the uncertainty in the distance does 
not dominate the uncertainties in the other parameters, we 



fix the normalisation of the neutron star atmosphere for a 
distance of 10 kpc. Spectral parameters determined by these 
fits are detailed in Table 0] The parameters we find for the 
first th ree Chandra observatio ns which were previously anal- 
ysed bv lWiinands et al] J2004) are consistent with the values 
these authors determined when assuming the same distance. 

2.2.4 Cooling curves 

From the spectral fitting we find that both and Fboi 
decrease significantly with time after the end of the out- 
burst. We measure an unabsorbed 0.5 - 10 keV flux of 
6.3I2 5 x 10~ 15 ergs cm -2 s _1 during the last Chandra ob- 
servation which is now lower than the upper limit from the 
previous non-detection with ROSAT'm the early 1990s when 
the source had been in quiescence for over 10 years. As 
with KS 1731—260 this cooling cannot be fit by a simple 
exponential decay curve, giving reduced x 2 values of 17.5 
and 13.6 for fits to the temperature and flux curves (see 
dotted line in Fig|^J . The cooling curves require an ex- 
ponential decay which levels off to a constant of the form 
y(t) = aexp [— (t — to)/b] + c. When fitting to the data, to 
was fixed to midday on the last day that the source was ob- 
served to be active, MJD 52159.5. The best-fitting cooling 
curves are shown in Fig. [3] For the curve a = 75.0 ± 2.8 
eV, b = 505±59 days, and c = 51.6±1.4 eV, with xl = 0.5. 
For the F bo i curve a = (4.5 ± 0.2) x 10" 13 ergs cm -2 s" 1 , 
b = 242 ± 13 days, and c = (1.7 ± 0.3) x 10~ 14 ergs cm -2 
s~\ with xl = 0.5. 

3 DISCUSSION 

We have presented Chandra and XMM-Newton observations 
monitoring two quasi-persistent neutron star X-ray tran- 
sients (KS 1731-260 and MXB 1659-29) in a quiescent 
state. The long accretion episodes onto both of these neutron 
stars before they went into quiescence significantly heated 
the neutron star crust out of equilibrium with the core. The 
monitoring observations cover a period of ~4 years since 
the sources returned to quiescence. Spectral fitting with 
an absorbed neutron star hydrogen atmosphere model to 
the data has clearly shown both of them cooling exponen- 
tially over time. From the in itial results of this monitoring 
JWiinands et alJl2002all2004l) and the preliminary results of 
IWiinandsl J2005I) it was not clear whether either of the neu- 
tron star crusts had reached thermal equilibrium with the 
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Table 4. Neutron star atmosphere model fits to the X-ray spectrum of MXB 1659—29 for 5 Chandra (CXO) and 1 XMM-Newton 
(XMM) observations. Luminosity is calculated assuming a distance to the source of 10 kpc. l-a errors on the parameters are given. The 
MJD given correspond to the mid-point of the observation. Note the column density parameter is tied between the observations when 
simultaneously fitting the data, and we determine Ng = (0.20 ± 0.02) X 10 22 cm -2 . 



ObsID MJD kT^g Bolometric flux Luminosity 

(Telescope) (eV) (10~ 14 ergs cm" 2 s" 1 ) (10 32 ergs s" 1 ) 



2688 (CXO) 


52197.8 


121 ± 2 


41 ± 3 


49 ±3 


3794 (CXO) 


52563.2 


85 ± 2 


10 ±1 


12 ±1 


0153190101 (XMM) 


52712.2 


77 ±1 


6 6+ - 4 
D - D -0.5 


7 q+0.5 

'■ 3 -0.6 


3795 (CXO) 


52768.9 


73+ 2 


r: 1+0.6 
°- i -0.4 


6 l +0 - 7 

"■'■-0.4 


5469 (CXO) 


53560.0 


58±l 


2 +0 - 4 


2 4+0.5 

z '^-0.6 


6337 (CXO) 


53576.7 


54+_i 


1.5 ± 0.5 


1.8 ±0.6 




1000 



Time since t Q (days) 



1500 



Figure 3. Cooling curves for MXB1659— 29. Top: Bolometric flux 
versus time since to , the last detection of the source accreting. The 
solid line shows the best fitting exponential decay to a constant. 
The constant is indicated as a dashed line. The dotted line shows 
the best fitting simple exponential decay, which does not fit the 
data well. Bottom: Effective temperature for an observer at in- 
finity versus time since the end of the outburst. The lines are as 
above. 



cores. However, the results of this latest monitoring indicates 
that KS 1731-260 and MXB 1659-29 have now reached 
equilibrium again and we have been able to measure the base 
level of the flux and effective temperature of these sources 
resulting from the state of the hot core. The additional data 
in this current work provides a significant improvement on 
the previous work with respect to constraints on the cooling 
curve. 



KS 1731—260 initially decreased by a factor of 2 in bolo- 
metric flux over the first half a year. Over the ~4 years since 
the source went into quiescence it has decreased in flux by a 
factor of ~5. It reached this level ~2.5 years after the end of 
the outburst and has remained at it since then, indicating 
that the crust returned to equilibrium with the core in this 
amount of time. This base level is set by the temperature 
of the core, which depends on the time-averaged mass ac- 
cretion rate onto the source over thousands of years. Fitting 
an exponential decay with a constant offset to the data we 
get the equilibrium temperature to be 70.0 ± 1.6 eV, and 
the base bolometric flux level due to thermal emission from 
the core as (9.2 ± 0.9) x 10 -14 ergs cm -2 s" 1 , which corre- 
sponds to a luminosity of (5.4 ± 0.5) x 10 32 ergs s" 1 . The 
e-folding times of the exponential decay were determined to 
be 325 ± 101 days for the temperature cooling curve and 
212 ± 60 days for the flux cooling curve. 

MXB 1659—29 decreased by a factor of 4 in bolometric 
flux over the first year after going into quie scence, and by a 
factor of 24 over the ~4 years of monitoring. IWiinands et alJ 
( 2004) found no evidence after the first 3 Chandra observa- 
tions of this source that the temperature or flux was reach- 
ing a equilibrium level set by the temperature of the core. 
The additional last 2 Chandra observations presented here 
indicate that the rate of cooling has decreased. Fitting an 
exponential decay with a constant offset to the data reveals 
that it has reached the base level which we determine to be 
51.6 ± 1.4 eV for the temperature and (1.7 ± 0.3) x 10~ 14 



ergs cm 



for the bolometric flux, which corresponds to 

|32 „„„„ „-l 



a luminosity of (2.0 ± 0.4) x 10 ergs s (assuming d=10 
kpc). The e-folding times of the exponential decay were de- 
termined to be 505 ± 59 days for the temperature cooling 
curve and 242 ± 13 days for the flux cooling curve. 

Rutlcdgc ct al ] J2002b) calculated detailed cooling 
curves for KS 1731—260 predicting the behaviour of the 
quiescent thermal emission from the neutron star given the 
source's accretion history (see their fig. 3). Their models 
compare how the emission could evolve for standard versus 
enhanced core cooling and low versus high crust thermal 
conductivity. These authors find that regardless of the core 
cooling, if the crust has a low thermal conductivity, then the 
crust should dominate emission over the core for ~30 years, 
after which time there is a large drop in the source luminos- 
ity. For high-conductivity crusts the timescale for this tran- 
sition is much shorter, ~1 year. Comparing these timescales 
to the observed cooling of KS 1731 — 260, then, in the con- 
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text of these models, the crust in this source must have a 
high thermal conductivity iWiinands et alj2002al : IWiinandsl 
2005), suggesting that it may be made of a pure iron lattice 
as increasing levels of impurity lower the conductivity. Un- 
fortunately no detailed cooling curves have been calculated 
for MXB 1659—29. As the cooling is dependent on the ac- 
cretion history of the source, we cannot directly compare 
the cooling of MXB 1659— 29 with the model curv es cal- 
culated for KS1731-260 by IRutledee et alJ ]2002bft . How- 
ever, in these models to see a rapid drop in flux within a 
few years of the source going into quiescence the crust must 
have a high conductivity, though there may be other phys- 
ical processes in the crust th at allow it to c o ol rapi dly. In 
particular, we note that the IRutledee et al ] J2002bft mod- 
els do not include Cooper-pair neutrino emission in the 
crust which will change the behaviou r significantly (e.g., 
Yakovlcv, Kaminke r fc Levenfisblll999ft . See the discussion 
in Ijonker et alJ (2006) for further comment on this differ- 



Once the crust has cooled, the quies cent luminosity 
is set by the emission from the hot core (IRutledee et alJ 
2002b). In the standard core cooling mechanisms the heat 
lost by neutrino emission is negligible and thus the heat 
deposited in the star is radiated as thermal emission. In 
the case where enhanced core cooling is allowed, the en- 
hanced level of neutrino emission, via the direct Urea 
process or Cooper-pairing, means that the vast major- 
ity o f heat deposited escapes as neutrinos (|C^h oL£LalJ 
l200ll: lUshomirskv fc Rutledeell200ll: IRutledee et alJ 12002b: 
iGusakov et al.l l2004). If enhanced core cooling occurs, then 
the quiescent luminosity set by the core temperature is 
significantly lower than if only standard core cooling oc- 
curs. We note that co ntinued cooling is predicted in the 
IRutledee et"ai] J2002b) models after the crust and core are 
in equilibrium. However, this cooling is at a significantly 
slower rate. 

Importantly, in both sources we have now measured the 
thermal equilibrium temperature, which is set by the core 
temperature, for the first time. The quiescent luminosity at 
the equilibrium level of both these sources is quite typical 
of the quiescent emiss ion from other 'nor mal' X-ray tran- 
sients (e.g. see fig. 5 in ljonker et al 1l2004al) whose outbursts 
are not as long as those of either source and last usually 
for only weeks to months. For the core of KS 1731 — 260 and 
MXB 1659—29 to have a similar temperature, and hence qui- 
escent luminosity, to the other X-ray transients then either 
the time-averaged mass accretion rate has to be similar to 
the other X-ray transients, or enhanced levels of core cooling 
are required. For the time-averaged mass accretion rates to 
be similar to other X-ray transients this would require ex- 
tremely long periods of quiescence of the order of hundreds 
to thousands of years if the observed outbursts are typical. 
For MXB 1659—29 there was only a period of 21 years be- 
tween outbursts and both known outbursts lasted for several 
years. The long-timescale accretion history for KS 1731—260 
is less clear with only one outburst having been observed, 
however, extremely long quiescent peri ods are hard to ex- 
plain with disk instability models (e.g. lLasotall200ll) . The 
measurement of the equilibrium lev el for both these sources 
confirms the previou s findings by JWiinands et alJ l2002al. 
120041 : Iwiinandsl2005ft from initial monitoring that enhanced 
core cooling processes are required for these sources. 



Comparing the e- folding timescales of the two sources, 
we find that KS 1731—260 cools at a faster rate than 
MXB 1659-29 (by a factor of ~1.5), suggesting that the 
heat conductivity of the crust in KS 1731 — 260 is greater 
than that in MXB 1659—29. This interesting difference is 
likely due to differences in th e properties of the cru st be- 
tween the two sources. In the (IRutledee et aljl2002bft mod- 
els this would mean that KS 1731 — 260 has less impurities 
in the crust compared to MXB 1659—29 as a pure iron crust 
has a higher conductivity than one containing heavier met- 
als also. However, one might expect the opposite to be the 
case, given that KS 1731—260 was actively accreting for 5 
times longer than MXB 1659-29. 

Given that we now have accurate cooling curves for 
these sources, and know the equilibrium level temperature, 
which places limits on the long-term time-averaged mass ac- 
cretion rate, theoretical models need to be fit to this data 
to determine how high the thermal conductivity needs to 
be and what levels of neutrino emissions are really needed, 
or any alternative properties that can account for the rapid 
cooling of the crust. As more than one outburst has been 
observed for MXB 1659—29, and because the cooling curve 
is much better constrained, we suggest that this source is 
possibly a better object to determine accurate theoretical 
cooling curves for than KS 1731—260. 

Although some neutron-star transients do not require 
enhanced core cooling to explain th eir quiescent lum i- 
nosity e.g. Aql X-l, 4U 1608-522 llColpi et alJ 1200 lTl . 
XTE J212 3— 058 (assuming th e recurrence timescale is 
>70 years, iTomsick et"ai] l2004h an d SAX J 1748. 8- 202 1 
in the globular cluster NGC 6440 dCackett et al.l 12005ft . 
there are several neutron-star transients in addition to 
KS 1731-260 and MXB 1659-29 that require either 
very-low time-averaged mass accretion rates or enhanced 
core cooli ng to explain th eir quiescent luminosity e.g. 
Cen X-4 dColpi et alJ l200ll). X1732-304 i n the globular 
clust er Terzan 1 Iw^nands et alJ 2002bJ^_JCacket^^dJ 
12006ft . RX J1709 30.2-263927 Ijonker et alJ l2003h . 
SAX J1810. 8- 2 609 jjonker et al.ll2004bl) , E XOl 747-214 
jTomsick et al] I2005T) and 1H 1905+000 (Ijonker et al] 
l2006h Whiiethe quiescent luminosity of KS 1731-260 
and MXB 1659—29 is comparable to many quiescent 
neutron star transients, it is significantly higher than 
the emerging group of quiescent neutron star transients 
that are very faint in quiescence, with quiescent lu- 
minosities of less than 10 32 ergs s . Although these 
faint quiescent sources include the millisecond X-ra y 
pulsars SAX J1808.4-3658 JCampana et al] 12002ft 
XTE J0929-314, XTE J1751-305 iWiinands et*ai]|2005bft 
and XTE J1807-294 JCampana et alJ 12005ft which may 
have different properties due to higher magnetic fields, 
there are also a number of 'normal' neutron-star X-ray 
transients that fall into this category (SAX J1810.8-2609, 
XTE J2123-058, EXO 1747-214). A notable faint quies- 
cent source is the quasi-persistent neutron-star transient 
1H 1905+000 as although it was seen to be in outburst 
for at least 11 years in the 1970s/80s, a recent Chandra 
observation did not detect the source in quiescence, with 

an upper limit of 1.8 x 10 31 ergs s " 1 (0.5-10 keV; for 

a distance of 10 kpc, Ijonker et al] 12006ft . Apart from 
KS 1731-260, MXB 1659-29 and 1H 1905+000, the only 
other quasi-persistent source to have gone into quiescence 
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is X1732— 304 in the globular clus ter Terzan 1. Using r ecent 
Chandra observations of Terzan 1 ICackett et alJ (120061) find 
that the most likely quiescent counterpart to this source 
has a 0.5-10 keV luminosity of 2.6 x 10 32 ergs s _1 , about 
6 years after the end of an outburst that lasted at least 
12 years. For comparison, the 0.5-10 keV luminosities of 
the last observations of KS 1731-260 and MXB 1659-29 
are 2.8lg;| x 10 32 ergs s _1 and 7St\% x 10 31 ergs s _1 
respectively. Therefore, KS 1731-260, MXB 1659-29 and 
XI 732— 304 all have higher quiescent luminosities than 
1H 1905+000. However, 1H 1905+000 was observed much 
longer after the end of its outburst then any of the other 3 
sources. While the difference in quiescent luminosity may 
suggest that these sources will continue to cool further 
over the next several decades, it is more likely due to 
differing long-term time-averaged mass accretion rates or 
neutron star properties. The differences in the properties 
of quiescent neutron star transients need to be explored 
further by additional observations to determine why some 
sources require enhanced core cooling and why some sources 
are much fainter than others. Observations of the next 
quasi-persistent neutron-star transient sources to return 
to quiescence (EXO 0748-676 and GS 1826-238 are both 
quasi-persistent sources that are still active) will help 
determine whether a decline in flux and temperature to 
a constant level on the timescale of years, rather than 
decades, is typical. 

We used an absorbed neutron star atmosphere model to 
fit to the data. If we use simple blackbody model in the spec- 
tral fitting instead, we get similar cooling curves, yet with 
different temperatures due to the differences between the 
shapes of blackbody and neutron star atmosphere spectra. 
When performing the spectral fits to both sources the nor- 
malisation of the neutron star atmosphere component was 
fixed to the distance of the sources. Unfortunately, these dis- 
tances are quite uncertain, thus, if the distances are wildly 
different from the values that we have assumed, the temper- 
atures and fluxes that we get from the spectral fits will be 
incorrect. However, if the normalisation is left free when fit- 
ting, then the uncertainty in this parameter dominates the 
errors in the other parameters. Importantly, the timescale 
of the decrease in the temperature and flux of these sources 
should be unchanged by a different distance, only the abso- 
lute values of the temperature and flux will change. If the 
sources are much closer than assumed then the temperature 
measured will be lower and the column density higher, and 
th e opposite if the source s are further away, as demonstrated 
bv lWiinands et all J2004) who fit the first three Chandra ob- 
servations of MXB 1659—29 assuming 3 different distances 
to the source ranging from 5-13 kpc. 

While we have interpreted the observations of 
KS 1731-260 and MXB 1659-29 in terms of the 
crust and core cooling model there are alternative 
models to explain quiescent neutron star emission 
including r esidual accretion or p ulsar shock emis- 
sion (e.g.. ICam pana et alJ Il9 98at iMenou et all Il999t 
ICampana fe Stellal2000HMenou fc McClintockl200ll) . While 
these are often used to explain the power-law tail present 
in the spectra of some neutron star transients, residual ac- 
cretion can produce soft spectra if accretion o ccurs onto 
the neutron star surface iZampieri et al] 1X9951. However, 
while variations in the luminosity of the sources in quies- 



cence can be explained by changes in the accretion rate or 
the amount of matter interacting with the magnetic field, 
our observations show that this would have to decrease 
exponentially in a smooth manner over a period of sev- 
eral years. Given that at the end of the outbursts of 'nor- 
mal' short-duration neutron-star tr ansients they reach their 
quiescent states w ithin weeks (e.g. ICampana et alj ri998b; 
Ijonker et" al. 2003) and that variations in accretion rate are 
generally not smooth, these alternative models seem un- 
likely. 

We note that from the initial observations of both 
sources it was found that a power-law component (which 
is often attributed to residual accretion) in addition to the 
neutron star a tmosphere component was no t required t o 
fit the spectra JWiinands et alJl2001al l2002al . l200l l2004h . 
While it has been observed that the contribution from 
the power-law component compared to the thermal com- 
ponent increases wit h decreasing quiescent luminosity be- 
low ~10 33 ergs s" 1 Jjonker et al.ll2004blal) . the statistics 
of the last few observations are unfortunately not good 
enough to check for this. However, if residual accretion 
is occurring during quiescence it should be on-top of the 
'roc k bottom' thermal emission due to deep crustal heat- 
ing ferown et al.lll998l : Etiitledge et alJl2002bl) that we have 
now detected. Several neutron star transients have been ob- 
served to vary during quiescence e.g. Aql X— 1, Cen X— 4 
and SAX J1748 8-202 1 i n the g l obular cluster NGC 6440 
jRutledee et alJ l2001allbl l2002al: ICampana fc StelTal 120031 : 
ICampana et al J 120041: ICackett et al]l2005l) . This variability 
can be interpreted as a variable power-law component, pos- 
sibly due to changes in the residual accretion rate or the 
amount of matter interacting with the magnetic field. Unfor- 
tunately this idea is hard to investigate with this data due to 
the low number of counts in the latter spectra. However, any 
variability in the flux from KS 1731-260 and MXB 1659-29 
over the next few years around the base set by the thermal 
emission from the neutron star may indicate that some resid- 
ual accretion is occurring. 
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APPENDIX A: RELATIONSHIP BETWEEN 
FLUX AND TEMPERATURE COOLING 
TIMESCALES 

Assuming that the X-rays are due to some blackbody-type 
emission, then the flux, F, and temperature, T, and just 
simply related via 

F = AT 4 , (Al) 

where A — j^a and R is the emitting radius, D the dis- 
tance, and a the Stefan-Boltzmann constant. We are mod- 
elling the cooling of the neutron star surface as an exponen- 
tial decay to a constant level, i.e., 

T(t) = Too + T'(t) (A2) 

where Too is the temperature at t = oo which we assume 
here to be constant, and T'(t) describes the evolution of the 
temperature of the neutron star crust with time, which we 
model as: 

T'(t) = aexp[-(t-t )/b], (A3) 

with a a normalisation constant, b the e-folding time, and 
to the start time. Combining Eg. lAll and IA2I leads to 

F(t) = AT(t) 4 = AlTco + T'(t)] 4 . (A4) 

At the start of the cooling, when Too <JC T'(t), then the 
T'(t) 4 term dominates, and Eg . IA4I becomes 

F(t) « AT'{t) 4 = Aae- 4(t - ta)/h . (A5) 

and there should be a factor of 4 difference between the e- 
folding times of F and T. This was f ound to be the case fo r 
the initial cooling of MXB 1659-29 iWiinands et alJl2004l) . 
However, further on in the cooling, once T'(t) ~ Too, then 
the e-folding time of F becomes a complicated function of 
the e-folding time of T due to the constant term. Finally, 
once the cooling is negligable and T'(t) <^ Too, then obvi- 
ously, F(t) = ATto- 
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